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Abstract 

The adiabatic cristal model is offered. It is shown that springy nuclei oscillations 
relatively electronic envelops and waves of such oscillations (inherent oscillations and 
waves) may exist in crystals. The analysis of experimental temperature dependencies 
of resistivity in semiconductors with electron- vibrational centres has shown that inherent 
oscillations effectively interact with crystalline phonons as well as with electrons and holes, 
creating powerful interaction of electrons and holes with phonons. The experimental 
narrow peaks of phonon-drag thermoelectric power at Debye temperatures from 77K to 
700K confirm existence of inherent oscillations waves in crystals. Inherent oscillations 
and waves gives rise to strong electron-phonon interaction and probably can bring about 
superconductivity at temperatures as below so and well above room temperature. 

Oscillations of cristals traditionally have always researched by analysing motion equa- 
tions of their models. Atoms in such models are replaced by spots with masses of atoms. 
These models do not correspond to an adiabatic theory of solids because the nucleus and elec- 
tron envelop of each atom in these models is considered to be hard are bound. Becouse in 
such conditions exchange of the energy between nuclei and electronic envelops is possible, the 
adiabatic principle is broken. That is why traditional crystal models are nonadiabatic. It's 
possible to find out new phisical properties of crystals by researching their adiabatic models. 
So substantiated adiabatic crystal model is important to build. 

1. Description of adiabatic model 

It's possible to build an adiabatic model of cristal in accordance whith adiabatic solution 
of stationary Shredinger equation: (T e + T z + V)ip = Wip, where T e and T z - kinetic energy 
operators of electrons and nuclei, V - a cristalline potential, ip - wave function and W - energy of 
cristal. This equation may be re-arranged to two equations supposing i/j = 0$ and separating 
variables: 

(T e + V)<j> = Ec/>, (1) 

(T 2 + E + A)$ = W<£>, (2) 

where E - energy of electrons and A - an adiabatic potential. Wave function describes electrons 
motion in the cristalline potential field V. Wave function $ describes nuclei motion. Because 
of presence adiabatic potential A in the equation (0) the functions and $ are mutually 
dependent, therefore electrons and nuclei in the cristal have a possibility to exchange the 
energy one with another. So the problem of cristal oscillations study on bases of equations ([I]) 
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and (0) in general is nonadiabatic. However, if potential A is small and its contribution to the 
energy of cristal can be neglected, it follows the adiabatic approach of Born-Oppenheimer [|I|]. 
In this approach an exchange of the energy between electrons and nuclei in cristal is absent 
and equation (|j) has a solution regadless of equation ([[]). Hartree-Fock method p], || may 
be applied to the equation (H) to define an efficient potential V(Rj), which is dependent on 
coorinates only of j-th nucleus, and reduce a problem about nuclei motion to one-partical 
problem: [Tj + V(Rj)]&j = Wj&j, where $j - wave function of nuclei, Tj - a kinetic energy 
operator, Wj - a spectrum of j-th nucleus stationary energy levels in the potential field V(Rj). 
Potential V(Rj) is defined by all electrons and nuclei exept j-th nucleus. Analysis shows that 
the main contribution to V(Rj) is introduced by s-electrons of K, L and M electronic orbits of 
j-th atoms. Minimum V(Rj) defines a position comparatively which nucleus is capable realize 
oscillatory motion. It is shown in Q that adiabatic method equitable if the energy of nucleus 
oscillations does not exceed energy of electronic transitions. Under such condition the energy 
of nucleus oscillations can not be sent to electrons. So electronic envelops stay nondeformed 
and moutionless but oscillatory motions of nuclei are principally possible. These oscillations 
can be called inherent (I) oscillations ||, since their properties are defined by inherent atom 
parameters - by mass, charge of nucleus and potential V(Rj) near centre of electronic envelop. 
Consequently, every atom of crystal in adiabatic model must be presented in a manner of 
inherent (I) oscillator which consist of nucleus and electronic envelop bounded with each other 
by springy force. 

Nuclei oscillations relatively atom envelops occur in small interval (~ 10~ 2 A°) and in general 
it is nesessary to research them by quantum methods. But, if study harmonic oscillations then 
it is sutable to use the known correspondence between results of quantum and classical theo- 
ries. This correspondence consists in coinciding the frequency of transition between quantum 
adjasent oscillatory levels of harmonic oscillator with his classical oscillations frequency. So 
energy spectrum of cristal harmonic oscillations is possible to reseach by clasical method. We 
used this possibility for description of oscillations in adiabatic simple univariate cristal model. 
Such model is shown on insertion fig. 1 where circles are electronic envelops of atoms and spots 
in circles centres are nuclei . Lattice spasing is equal to a. System of motion equations for this 
model may be written as: 

m^K = -$A<-K)\ (3) 
= -eiK - K) - - K-x - K + il (4) 

where M - nucleus mass, m - an electronic envelop mass, U' and U" - displacings of nucleus 
and electronic envelops, Q\ and 82 - force factors, t - time and elementary cell number n = 
0, ±1, ±2, ■ ■ • . If search solutions of this system of equations as harmonic waves, so the 
dependency of angular frequency u from wave vector q can be expressed as 

wi, 2 (q) = (G/2){1 ± (4:9i9 2 )/MmG 2 )}, (5) 

G = (3/m* + •yC/m, C = 4sin 2 (aq/2), a - lattice spasings, m* = (1/M + l/m) -1 . Qualitative 
curves u;(q) are shown on the fig. 1 and contain known acoustic (A), as well as new (I) branches. 
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In complex real cristals except known acoustic and optical branches there are also branches of 
inherent oscillations. The number of I-branches is equal to 3r (r - a number of atoms in the 
utit cell of crystal). I-oscillations and waves can exist even though elementary cell kepts only 
one atom and optical oscillations are absent. 



2. Distinctive properties of iherent oscillations and waves 

Energy spectrum of inherent oscillations is possible to calculate taking into account an 
interaction between oscillations of different electronic envelops. Displacements of electronic 
envelops may be coherent within certain area with the distinctive size A due to the interaction 
between them. Thermodynamic analysis of adiabatic crystal model gives the value: A = 



Zm*/ne 2 /i, where Z - an atomic number, m* - effective mass of electronic envelops, e - a 
charge of electron, n - electronic density, /i - magnetic constant ||. The value A may exeed 
groups of ten lattice spasings. In such conditions a nucleus of each atom can realize oscillatory 
motion relatively big mass of crystal coherent area. Analisis shows that a, (3, 7 - types of 
inherent oscillations and waves are to exist in dependence of displacements K and L - orbits, 
ct-tipe of I-oscillations corresponds to nucleus oscillations relatively electronic envelop. Joint 
nucleus and K - orbit oscillations relatively remaining part of electronic envelop correspond 
to (3 - oscillations. Joint nucleus, K and L - orbits oscillations relatively remaining part of 
electronic envelop correspond to 7 - oscillations. Elementary quantum (huj z ) of a - oscillations 
for neutral atom with number Z > 2 calculated with provision for shelding of nucleus by 
electrons can be expressed as: 



Tiuj 7 



e 2 Z-(-n 3 X (Z-0 



(6) 



\ 3ire (m n + m p ) \ \a J Z-( j Z 

Shelding is taken into account by ( = 5/16 and £ = sZ 1 ^ 3 , s ranges from 1 to 1.15 when changing 
an atomic number Z from 8 to 80, x = 1-2 takes into account a contribution in electronic density 
from s-states of L, M and N - orbitals, $ = 0.88534, ao = 0.529A , Eq - electrical constant, m n 
and m p - masses of neutron and proton. Elementary quantum of (3 - oscillations possible to 
define on the eq. (|6|), having place x — 0.2. Elementary quantum of 7 - oscillations possible 
to define on eg. having place x = 0.056. Energy spectrum of harmonic inherent a, (3 or 
7 - type oscillations and waves are possible to describe by harmonic oscillator formula 

= + 1/2), (7) 

where oscillatory quantum number v = 0, 1, 2, • • •, and E = huj z for oscillations a, (3 or 7-type. 



3. Interaction of inherent oscillations with phonons 

Inherent oscillations and waves may exist in the ideal crystal but generate such oscillations 
and waves is possible (to the account of electrons or holes energy) by means of local centres with 
strong electron-phonon interaction. Inherent oscillations and waves distort a crystal. They are 
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capable effectively interact with other springy oscillations as well as with electrons and holes, 
ensuring efficient electron-phonon coupling. This can cause a phonon-drag component of ther- 
moelectric power and other physical effects. So interaction of inherent oscillations with acoustic 
phonons is possible to track at simple univariate crystal model adding in to the equation (|4]) 
additional force -SU£. This force describes effect from the coherent areas. Corresponding de- 
pendencies o;(q) is possible to get by replacing G on (G + 5U") in eq. ([|). Acoustic branch 
suffers most alteration alongside centre of Brillouin area. At the condition 5 > a forbidden 
frequency area ■ • • uj* appears. At condition 5 < a forbidden area for wave vectors ■ ■ ■ q* 
appears. Dependencies cu(q) for 5 > and 5 < are shown qualitatevely on fig. 1 by dashed 
qurves. Obviously at the conditions uj* > and \q*\ > experimental acoustic methods for 
crystal oscillations studies may turn out to be ineffective. 

Acoustic phonons at the conditions < uj < uj* or < |q| < |q*| can not exist in crystals 
and consequently can not scatter charge carriers. For this reason superconductivity is possible 
if other mechanisms of scattering are absent. 

4. Experiments, results, discussion 

It is known that at temperatures T < 70K thermoelectric power (TEP) may contain a 
phonon-drag component (PDC). Absence of PDC at T > 70K earlier was explained by in- 
sufficiently strong coupling between electrons and phonons ||. Such opinion in common was 
saved up to present time. However TEP in cristalline ropes of carbon nanotubes supposedly 
was explained by contribution of PDC at temperatures from 4.2K to 300K [^]. In carbon 
nanotube films on substrates at temperatures up to 600K narrow pius of TEP are discovered, 
connected with PDC || |J. PDC in carbon nanotube films on substrates is stimulated by 
the interaction of electrons in film with phonons of substrate due to inherent a-oscillations of 
Carbon (E = 0.25 eV) and Oxygen (E = 0.22 eV) atoms. In this work we measured PDC in 
monocristals GaP and Si, in porous silicon on silicon substrates. 

In experiments were used flat semiconductor samples with the thickness 200 mkm contain- 
ing local centres with the strong electron-phonon coupling. We investigated GaP samples with 
impurities of aluminum or sulphur 10 15 cm -3 ): GaP(Al) and GaP(S). Such impurities were 
choosen becouse atoms Al and S have masses vastly exeed mass of atom Ga. This advantages 
emergence electron-vibrational centres and generation of inherent oscillations. Investigated sili- 
con samples had impurities of Phosphorus (~ 5 • 10 15 cm -3 ) and Oxygen^ 10 18 cm~ 3 ): Si(P,0). 
Before measurements the samples where subjected to heating in vacuum within 5 minutes 
at T=600K then they were cooled up to room temperature within 0.2 minues. Concentration 
of Oxygen in Si was determined on the basis of data about intensities of distinctive optical 
absorbtion near 9 mkm. Layers of porous silicon (Si*) with the thickness about 0.3 mkm on 
n-silicon substrates with resistivity 3 Ohm- cm were also investigated. In samples Si(P,0) and 
in Si* local centres seem to be formed by atoms of oxygen (A-centres) for which the constant 
of electron-phonon interaction is equal to 5. 

Their were measured temperature dependencies of electrical resistivity p(T) and tempera- 
ture dependencies of thermoelectric power in samples within the range from 77K to 700K for 
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the reason of studying the influence of inherent oscillations and waves on migration of charge 
carriers. Temperature difference at mesurements of TEP did not exeed 3K ± 0.2K. The spec- 
trums of infra-red reflection modification (dR) were measured, caused by impurities in GaP 
samples, in the optical range from 15 mkm (83 meV) to 2 mkm (620 meV) at 300K. The angle 
between nonpolarized light beam and flat surface of sample was equal to 45°. 

Typical temperature dependencies of electrical resistivity p(T)/p of GaP samples are 
brought on fig. 2 in semilogarithmic coordinates, po is a constant and does not influence 
magnitudes E a . Magnitudes p are selected for each curve that the curves were conveniently 
plased on the fig. 2. Cuve 1 on fig 2 referes to temperature dependencies of electrical resistivity 
of GaP sample without impurities. Typical dependencies p(T)/p of GaP(Al) and GaP(S) sam- 
ples are presented on fig. 2 by cuves 2 and 3. These curves are broken straight lines. Tangent 
to these curves select areas which correspond to concrete values E a . Values E a , measured at 
T < 330-fT, are contributed in Table 1 where stars choose the lines with energies fit known ener- 
gies of phonons in GaP [|10j which are brought here in the central column. Activation energies 
of GaP(Al) and GaP(S) measured at T > 330/T are contributed in Table 2 where values of 
oscillatory energy of impurity atoms olso brought, calculated on the eq. (0) with provision for 
values (E = 0.283eV) for Al and (E = 0.301eF) for S. Modification of optical reflection (dR) 
caused by impurity atoms Al in GaP are presented on fig. 3. Experimental temperature depen- 
dencies of TEP for GaP samples are presented on fig. 4. Temperature dependencies of TEP for 
samples Si(P,0) and for polisilicon film on the silicon substrate (Si*) are presented on fig 5. 

Cuve 1 on fig. 2 may be characterized with activation energy E a = 0.7 meV at temperatures 
T<330K but value E a is close to the width of GaP forbidden gape (2.4 eV) at temperatures 
T > 330i^. Curve 2 and 3 on fig. 2 may be characterized by several activation energies. 
Corresponding values E a are contained in Table 1 and Table 2 and may be explained differently. 
Values E a in lines of Table 1, marked by the stars in the first column, close to crystalline phonons 
which intensively interact with local centers. These phonons are pointed out in central column 
in Table 1. It is impossible to explain the experimental dependencies p(T) by dissipation 
of charge carriers by phonons since the dissipation is capable to create an opposite effect to 
the observed effect of reducing resistance when increasing the temperature. We connect these 
energies with generation of free charge carriers from electron- vibrational levels of local centers 
created by atoms Al or S. Other energies in Table 1 supposedly may be explained by generation 
of free charge carriers from electron-vibrational levels formed by inherent oscillations (j3 and 
7-types) of Al or S and possible by combinations of such oscillations with cristalline phonons. 

It is seen from Table 2 that activation energies of samples with each type of impurities can 
be divided into two groups which corresponds to two columns dextral. One group consists of 
activation energies which are described by the equation for quantum harmonic oscillator (|7]). 
These activation energies are connected with inherent oscillations of impurity atoms and cor- 
respond to transitions from oscillatory energy levels with z/ = 0, 1, 2, •••in the minimum 
of oscillatory potential where oscillatory energy is zero. Such transition for free quvantum 
harmonic oscillator are forbidden but they are possible for classical oscillator. Consequently, 
inherent oscillators of impurity atoms show duality of properties that can be explained by their 
interaction with the cristal. Other group of energies in Table 2 consists of energies multiple 
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Table 1: Activation energies for GaP(Al) and GaP(S) at T < 330i^ 



Activation energies E a (meV) for 
GaP(Al) samples with numbers: 


Phonons in 
GaP (meV) 
PI 

II II 


Activation energies E a (meV) for 
GaP(S) samples with numbers: 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


8.1 






7.2 


8.0 














* 15.0 






14.8 


14.5 


14.25 


15.0 


14.8 






14.3 


* 24.5 








25.0 


24.42 


24.0 




24.6 








35.0 
















28.0 


28.0 


* 45.0 








44.6 


44.75 




42.0 








* 




49.0 


48.0 




47.00 


47.0 










75.0 




83.0 












70.0 







Eq. Energies of this group are also connected with inherent oscillations of impurity atoms 
and correspond to transitions between different oscillatory energy levels (between levels with 
different v). It is seen that the value Eq is the same for both groups of energies. Consequently, 
both groups of energies refer to the same type of centres that show a classical and quantum 
properties (duality of properties). 

Studies of infra-red reflection (R) also confirm presence of inherent oscillations in GaP(Al) 
and GaP(S) samples and their intensive interaction with electrons and crystalline phonons. 
Modification of infra-red reflection spectrum (dR) caused by impurities atoms was separated 
on components in accodance with complex oscillatory model [11, 12] taking into account the 
contributions from j different charged oscillators. Theese components are marked on fig. 3 by 
numerals j: 2, 3, 4 and 5. Every component involved contribution into reflection (R) reach the 
highest value when optical frequency (u) satisfy to the conditions u p > uj > Q. Q - an oscillator 
frequency and io v - a frequency of springy oscillations of lattice. Minimum dR locates near uj p . 
Agreement of experimental - 1 and amount of calculated spectrums 2, 3, 4 and 5 is reached 
if energies Tiui p are equal to the energies of a - tipe inherent oscillations for Al: 0.5Eq, Eq, 
1.5-Eoj 2i?o- Two from them comply with calculated on eq. (0) if Eq = 0.283 eV and v = or 1 
but two other are multiple the same E = 0.283 eV. Reflection spectrums of GaP(S) also are 
described within the complex osacillator model jn], [T2|] if Eq = 0.301 eV. Energies HQ for both 
types of impurities (A1,S) may be connected with 7 - type of inherent oscillations Al(61.1 meV) 
or S(65.0 meV). Energy loss of inherent oscillators is great (G/fl ~0.09; G- damping factor) 
and corresponds to the tight binding of impurity atoms with the cristalline phonons. Thereby 
centres formed by impurities Al and S in GaP show duality of optical properties. Theese 
properties are defined by interaction of inherent oscillations of impurities atoms, crystalline 
phonons and electrons with each other. 

The best consent of calculated and experimental reflection spectrums is reached if optical 
dielectric permeability (e = 2) is small in contrast with the highfrequency dielectric permeability 
GaP (e = 8.457) |fl3| . Probably measured value of dielectric permeability follow refere to local 
centres in which electronic optical transitions occur. A cross section of photons capture can be 
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Table 2: Activation energies for GaP(Al) and GaP(S) at T > 330K 



E a (eV) for GaP(Al) 
samples whith numbers: 


Calculated 
on eq. (171) 

1 \l 1/ 


Multiple E 

(eV) 


1 


2 


3 


4 


5 






0.14 


0.14 


0.14 


0.138 


0.137 


0.1415 = 0) 




0.28 


0.29 


0.29 


0.28 


0.28 




E = 0.283 


0.42 


0.42 


0.43 


0.42 


0.43 


0.4245 (u = 1) 






0.57 




0.56 


0.58 




2E = 0.566 


0.71 




0.72 






0.7075 (u = 2) 










0.85 






3E = 0.849 




0.97 








0.9905 (u = 3) 




1.10 








1.11 




AEq = 1.132 


E a (eV) for GaP(S) 
samples whith numbers: 






6 


7 


8 


9 


10 






0.15 


0.15 


0.15 


0.15 


0.15 


0.1505 [y = 0) 




0.3 


0.29 


0.3 


0.3 


0.31 




E = 0.301 












0.4515 {v = 1) 




0.6 




0.6 


0.61 






2E = 0.602 










0.74 


0.7525 (v = 2) 








0.92 








3E = 0.903 


1.03 










1.0500 (v = 3) 





defined by wavelength of cristalline phonons which interact with the centres. 

Inherent oscillations can spread in crystals as waves of inherent oscillations. Inherent oscil- 
lations, crystalline phonons and electrons (holes) powerfully interact with each other and give 
rise to electrical currents. This is confirmed by particularities of TEP. Curve 1 on fig. 4 is typical 
temperature dependency of TEP for samples GaP(S). Curve 2 on fig. 4 is typical for GaP sam- 
ples without impurities and may be explained by close to inherent conductivity. Temperature 
dependency of TEP for GaP(Al) is like curve 1. Curve 1 on fig. 4 keeps pius which marked by 
arrows and by letters. Polarity of these pius agrees with the polarity of drift TEP. We connect 
these pius with the phonon drag of electrons. Pius A, B, C and F are situated at Debye temper- 
atures of cristalline phonons in GaP Q: 95K (TA,8.2 meV); 168K (TA, 14.25 meV); 288K (LA, 
24.42 meV); 542K (LO, 44.75 meV). Broad pius D (» 345K) and E (» 475K) can to be ex- 
plained by combinations of cristalline phonons: (TA+TA, 28.6 meV) and (TA+LA, 38.67 meV). 
Investigations of p(T), TEP temperature dependency and infra-red reflection show that inher- 
ent oscillations Al and S in GaP actively interact with cristalline phonons and electrons. This 
ensures electron-phonon interaction sufficient for the realization of phonon drag effect (and 
probably superconductivity) at very high temperatures. 

Curve 1 on fig. 5 represents peculiar temperature dependencies of TEP for samples Si(P,0). 
It has pius spesified by arrows and by letters a, b, c, d, e. Polarity of these pius coincides 
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with the polarity of drift TEP. Pius a, b, c are situated at Debye temperatures of acous- 
tic phonons with wave vektors along certain directions [I3|: <111>, 200. 4K (16.7 meV); 



<110>, 214.8K (17.9 meV); <100>, 252K (21.0 meV). We explain these pius by phonon drag 
of electrons. Pius d and e we involve with TO phonons in Si. Curve 2 on fig. 5 represents the 
temperature dependency of TEP in Si* and has pius p, q, r, h with different polarities. Tempera- 
tures of these pius are equal to Debye temperatures of phonons in critical points in the Brillouin 
area of Si: L(W) - 551K (45.9 meV); L(L) - 606K (50.5 meV); TO(X) - 683K (56.9 meV) and 
TO(L) - 712K (60.9 meV) accodingly. We explain pius f, g, h, i by phonon drag of holes but 
peak j we explain by phonon drag of electrons. 

The sufficient interaction between electrons (houls) and phonons in silicon samples can 
be explained by inherent oscillations of centres which are formed by impurities of Oxygen 
(E =0.22 eV). These are to be A-centres. Such interpretation agrees with experimental tem- 
perature dependencies of electrical resistivity for the same samples. These dependencies are 
discribed by several activation energies among which there are multiple E$ = 0.22 eV or equal 
to calculated on the eq. [7] if E is equal to the same value. 

5. onclusion 



The analysis of experimental results has shown that phonon drag effect is caused by the 
phonons which actively interact with electron-vibrational centres formed by impurity atoms 
(Sulphur and Aluminium in GaP, Oxygen in silicon samples). Such centres can generate inher- 
ent oscillations and waves. Theoretical possibility of strong electron-phonon interaction, caused 
by inherent oscillations and waves, was confirmed by the phonon drag effect at temperatures 
up to about 700K. The phonon drag component of thermoelectric power also has confirmed the 
presence of inherent oscillations waves in crystals. The expected superconductivity from about 
room temperatures up to temperatures more 700K provided by inherent oscillations and waves 
probably will be realized. 
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Figure captions 

Fig. 1. Inherent (I) and acoustic (A) branches of oscillations in adiabatic simple univariate 
crystal model ( 1 - at 5=0; 2 - at 5 >0; 3 - at S <0 ). Fragment of adiabatic simple 
univariate crystal model is shown on insertion. 

Fig. 2. Temperature depenencies of electrical resistivity p(T) for samples GaP: 1 - without 
impurities, 2 - with impurity of Aluminium; 3 - with impurity of Sulphur. 

Fig. 3. Change of optical reflection (dR) caused by impurity of Aluminium in GaP — 1. 

Calculated components of dR connected with inherent oscillations of Alu- 
minium : 2, 3, 4, 5. 

Fig. 4. Temperature dependency of thermoelectric power: 1 - in GaP with impurity of 
Sulphur; 2 - in GaP withaut impurities. 

Fig. 5. Temperature dependency of thermoelectric power: 1 - in Si with impurity of Phos- 
phorus and Oxygen; 2 - in porous silicon (Si*). 
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